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1. Motivation 


In an isospin symmetric world, the up (n) and down {d) quarks are identical particles. It is 
known {cf. Table 1) than in Nature isospin symmetry is explicitly broken by the non-zero mass and 
electric charge differences of the u and d quarks. However, the effects of this breaking are expected 
to be small relative to typical strong interaction energies such as hadron masses. Indeed, it is clear 
that the light quark mass mass difference 5m = ntu — mj represents one percent or less of any typical 
QCD energy scale. Similarly, the strength of the electromagnetic (EM) interaction relatively to the 
strong one is essentially given at low energy by the fine structure constant a ~ 0.007. For those 
reasons we can reasonably state that, for observables with a non-vanishing isospin symmetric part, 
isospin symmetry is a good approximation of reality with an 0(1%) relative error. 

Nevertheless, it is interesting to note that these small isospin breaking corrections are crucial 
to describe the structure of atomic matter in the Universe. Indeed, one particular effect of isospin 
symmetry breaking is the mass splitting between the proton (p) and the neutron (n). This mass 
difference is known experimentally with an impressive accuracy [1]: 

AMa,=M„-Mp = 1.2933322(4) MeV (1.1) 


The sign of this splitting makes the proton and the hydrogen atom stable physical states. Also, the 
size of AMn determine the phase space volume for the neutron j8-decay n ^ p + e^ + Ve- At early 
times of the Universe (t ~ 1 s and T ~ 1 MeV) and under standard assumptions', the existence of 
-decay allows to infer that the ratio of the number of neutrons and protons is approximatively 
equal to: 


~ exp 


(A^) 


( 1 . 2 ) 


This ratio is one important initial conditions of Big Bang Nucleosynthesis. Also, in our actual 
Universe, j3 -decay and its inverse process are known to be responsible for the generation of a large 
majority of the stable nuclides chart though nuclear transmutation. Even if the nucleon isospin mass 
splitting is a well known quantity, predicting it from first principles is a difficult problem because 
of the complex non-perturbative interactions of quarks inside the nucleon. The proton carTies an 
additional EM self-energy compared to the neutron, so just from QED one would expect to have 
AMn < 0. However, the fact that the experimental value of AMn has the opposite sign suggests that 
the strong isospin breaking effects are competing against the EM effects with a larger magnitude. 
This would mean that an important part of the structure of nuclear matter as we know it relies on 



u 

d 

Mass (MeV) [1] 
Charge 

2.3(1°-) 

4.8(1°-) 


Table 1: Physical properties of the up and down quarks. 

*The neutrino number density ny/ny is assumed to have the order of the baryon density number which is very small. 
This assumption is not valid anymore in some new physics scenarios but even in these hypothetical cases rin/rip depends 
strongly on Mljy. 
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a subtle cancellation between the small EM and strong breaking effects of isospin symmetry in the 
nucleon system. Therefore, it is fundamental to have a theoretical understanding of the nucleon 
isospin mass splitting. 

Considering that isospin breaking effects in the hadron mass spectrum are generally measured 
quite precisely, it is also interesting to understand how one can use this information to deduce the 
masses of the individual u and d quark masses. For example, it is important to know if = 0 
could be a realistic solution to the strong CP problem. While recently (c/ the FLAG review [2]) 
considerable progress has been made in determining precisely the average up-down quark mass 
niud from first principles, such a computation is still missing for the individual masses. Because 
the kaon {K) is a pseudo-Goldstone boson of chiral symmetry breaking, the isospin mass splitting 
AM|- = M|.+ is very sensitive to 5m. But in order to extract 5m, one has to understand how 

to subtract the EM contribution to this splitting. One well known result in this direction is Dashen’s 
theorem [3] which states that, in the SU(3) chiral limit, the EM Kaon splitting is equal to the EM 
pion (tt) splitting: 

AgEoAfl- = AqedMI + 0{ams) (1.3) 

This result is important because it is known [2] that with good accuracy, AqeuM^ ~ AM^. The 
remaining question is: how large are the 0{ams) corrections in (1.3)? One way to quantify these 
corrections is to consider the dimensionless quantity £ defined in [2] as follows: 

_ AqedM|- - AqedM^ 

AMI ^ ^ 

This quantify is consfrucfed such fhaf if vanishes in the SU(3) chiral limit. There were several 
attempts in the 1990s to compute these corrections analytically from effective theories which leaded 
to controversial results. It makes this quantity a good target for a lattice calculation. 

The problems presented in this section, and more generally in any computation of isospin 
corrections to low-energy QCD observables, are difficult to solve because of the highly non- 
perturbative behavior of the strong interaction in this regime. It has been shown [2] that it is now 
possible to predict fundamental isospin symmetric QCD observables through lattice QCD simula¬ 
tions with a full control over the method’s uncertainties. It is then reasonable to think that lattice 
simulations could be a reliable way to understand and compute isospin breaking effects. Moreover, 
besides the physical interest of these effects, actual lattice calculations are reaching a sub-percent 
precision on several standard observables and the assumption of isospin symmetry is becoming the 
dominant source of systematic uncertainty. 


2. Lattice QCD+QED 

In this section we review how to add EM interactions to lattice simulations. As we will see, 
the main difficulties comes from the singular infrared structure of QED. We explain a possible way 
to define QED in a finife volume and whaf are fhe associafed finife-size effecfs. We fhen discuss 
fhe discrefizafion of fhe fheory and fhe simulation fechniques used so far. 
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2.1 QCD+QED in a finite volume 

Let us consider a diagrammatic contribution ^ to a correlation function featuring a photon 
loop {e.g. the 1-loop part of the electron EM self-energy). In infinite volume, will have the 
following form: 

i d^A: 1 

^oc = y (2.1) 

The integral (2.1) may be ultraviolet (UV) divergent, which can be dealt with through renormal¬ 
ization. The photon pole at = 0 can also generate infrared (IR) divergences. However, (2.1) is 
not mathematically undefined per se, the undefined A:^ = 0 value of fhe infegrand is jusf a poinf 
(sef of measure 0) and can be ignored. Moreover, if is known fhaf in some cases {e.g. fhe on-shell 
self-energy of a particle), fhis singularify is infegrable. In a finife volume wifh femporal exfenf T 
and spafial exfenf L, momenfa become quantized in a way depending on fhe choice of boundary 
conditions. If one chooses periodic boundary condifions for fhe phofon field, fhen fhe confribufion 
(2.1) becomes: 

Q{T,L) = ^ ^ ^f{k,pi,...,p„) (2.2) 

'' ^ keBZ{T,L) ^ 

where: 

BZ(r,L) = ^Zx^z3 (2.3) 

Now fhe expression (2.2) has an isolafed, undefined confribufion coming from fhe phofon pole 
which cannof be summed in any way. As menfioned in [4, 5], fhis singularify is classically relafed 
fo fhe facf fhaf Gauss’ law does nof aufhorize a nef charge fo exisf in a finife, periodic volume. 

2.1.1 Photon zero-mode subtraction schemes 

If one wants to keep periodic boundary conditions on the photon field, one possible solution 
to deal with the zero-mode singularity is to remove a subset of mode containing 0 from the finite- 
volume degrees of freedom. This will of course alter the physics in finite volume, but if the chosen 
subset converges in the infinite volume limit to a set of measure 0 then naively the physics in infinite 
volume remains unchanged. This is only a necessary condition, in principle one needs to check that 
the subtraction procedure does not accidentally couple the IR and UV structure of the theory which 
could introduce a complicated volume-dependent renormalization of the theory. 

Naively, the most minimal zero-mode subtraction procedure is to set the k = 0 mode of the 
photon field to 0. Following [6], we name the resulting theory QED-pL- This finite-volume pre¬ 
scription has been used for numerical calculations in [7, 8, 9]. Although this scheme is simple, it 
introduces some strong finite-volume effects which can be hard to control. Indeed, considering the 
T —)• -|-oo limit at fixed L of the QEDjl version of (2.2), one obtains: 

^qeDtl(T,T) = ^ ^ f{k,pu...,p„) ^ £ f{k,pi,...,pn) (2.4) 

'' ^ teBZ(T,L) ■' J Zn |5gg2{L) 

kfO 

where BZ(L) = '^1?. Because it is one-dimensional, the integral in (2.4) might be IR divergent 
even in cases where its four-dimensional version converges. As an example, finite-volume effects 
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on the 1-loop mass correction in spinor QEDjl were computed in [6]: 


mQED-^i^{T ,L) 


= m ■ 




f. 2 

\ ^ 

(. 2 


% r1 

< 1 -q^a 

2mL 


1 - 



Stt 

{mL)~ 


1 - 


coth(mr) 


371 L 
2{niLY T 


(2.5) 


where m is the infinite volume mass, q is the charge in units of e and K = 2.83729... is a known 
numerical constant. This expansion is exact up to corrections that decay exponentially in the in¬ 
finite volume limit. In this example we explicitly see a term proportional to ^ which represents 
the IR divergence related to the limit (2.4). In conclusion, it appears that QED-pL has two cumber¬ 
some properties. Eirstly, the infinite-volume limit has to be taken with special care (i.e. by keeping 
bounded) and secondly this extrapolation depends on the aspect ratio As discussed in [6], 
the singularity of QED^l in the large T limit can be explained in the following way. The photon 
zero-mode removal can be implement by adding a non-local term in the Eagrangian of the theory. 
This term couples values of the electromagnetic potential on different time-slices, breaking the 
reflection positivity of the action. So strictly speaking QED^l does not admit a quantum mechan¬ 
ical description and the divergence in T is a symptom of the lack of a thermodynamic limit. This 
singular behavior has been discovered independently by the MIEC collaboration [10]. 

An alternative to QED^l is to remove all spatial zero-modes, i.e. to set to zero all the photon 
modes k with k = 0. This scheme is inspired from [11] where QED is formulated in a finite 
spatial volume directly with an infinite temporal dimension. We denote this prescription QEDl [6]. 
Because it does not couple field values on different time-slices, QEDl has positive reflexivity and 
therefore a correct particle physics interpretation. In this theory, the finite-volume effects on the 
1-loop mass correction of a spin ^ particle are given by [6]: 


mQED, {T, L) = m\\-q^a 

Compared to (2.5), this relation is now completely independent from the aspect ratio QEDl has 
been used in electro-quenched simulations in [12] and in full simulations in [6]. 

The mass corrections in QEDl QEDjl were compared to numerical quenched QED 
(which is exact at the 1-loop) simulations in [6] and perfect agreement is found between the simu¬ 
lations, (2.5) and (2.6). These results are summarized in Eigure 1. 


K / 2 \ 371 

2mL I mL J {mL)^ 


( 2 . 6 ) 


2.1.2 Finite-volume effects on hadron masses 

All the previous statements were made for elementary particles which interact only through 
QED. Here we review how to generalize this discussion for hadrons, i.e. composite bound states of 
the strong interaction. All of the work presented here are computations performed in QEDl with 
an infinite time dimension. 

To study finite-volume effects on hadrons masses, one possible approach is to use low-energy 
effective theories of the strong interaction coupled to QED. This was done first for meson masses in 
the context of SU(3) partially-quenched chiral perturbation theory in [11]. The results of that work 
were studied numerically and generalized to an SU(2) plus heavy kaons theory in [12]. EM finite- 
volume corrections to meson, baryon, nuclei masses and to the hadronic vacuum polarization were 
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Figure 1: Finite-volume corrections to an elementary fermion mass in QED. The black (resp. red) points 
represent QEDxl (resp. QEDl) quenched QED simulations (which is exact at 0(a)). The black (resp. red) 
curves represent the theoretical prediction (2.5) (resp. (2.6)). The only fit parameter is the infinite-volume 
mass. The dashed blue line is the prediction from [5], the disagreement between this formula and the data is 
commented in section 2.1.2. 


studied in [5] using non-relativistic effective theories. In this context, the finite-volume corrections 
appear as the elementary particle ones plus terms depending on the structure of the particle (radius, 
polarizabilities, ...). 

However, there is a disagreement between the point-like limit of [5] and the QED predic¬ 
tion (2.6) concerning the fermion mass correction. The difference is a relative factor of 2 in the 
0[a/(/ML)^] correction. The QED simulations presented in Eigure 1 strongly favor (2.6). An ex¬ 
planation for this discrepancy has been recently proposed in [13]. In the non-relativistic limit, the 
particle and antiparticle degrees of freedom decouple and therefore in non-relativistic theories one 
does not expect the antiparticle modes to contribute to the particle mass corrections. However, as 
pointed out in [13], these modes contribute through the subtraction of the photon zero-modes in 
QEDl- Once properly added, this residual fermion-antifermion interaction generates a O [a/ {mL)^] 
finite-volume correction which solves exactly the discrepancy with (2.6). 

The structure of finite-volume corrections to hadron masses was also studied beyond the effec¬ 
tive theory level [6]. In that work it is confirmed fhaf fhe fwo firsf orders of fhe finife-volume expan¬ 
sion in ^ are universal and identical fo fhe pure QED case. These ferms are defermined by gauge 
invariance, fhi'ough fhe consfrainfs on fhe elecfromagnefic form facfors provided by fhe Ward- 
Takahashi idenfifies, and follow from fhe analyficify properfies of 1-particle iiTeducible Green’s 
functions in fhe relevanf quanfum field fheories. The sfrucfure confribufions only enfer af leasf af 
order 0[a/(mL)^] which can be seen explicifly for fhe effective fheories presenfed in [11, 12, 5]. 
The universalify of fhe fwo firsf orders is an importanf information as if allows fo impose fhese 
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corrections analytically in lattice data analyses without introducing any model dependence. 

2.2 Lattice QED 

There are essentially two approaches to discretize QED: a naive, non-compact action where 
the gauge potential is still the field variable or a compact Wilson action similar to lattice QCD. 
By naive discretization we mean that the lattice action is defined as follows: 

of 

~ ~Af ^ [dpAy[x) — dyAn[x)] (2.7) 

where a is fhe laffice spacing and is some firsl-order finife-difference operafor. So far, excepf 
fhe sfarfing projecf of fhe MILC collaboration [14], only fhe non-compacf action has been used in 
fhe confexf of laffice QCD-i-QED simulafions. One of fhe main mofivafions in making fhis choice 
comes from fhe facf fhaf fhe non-compacf action is free and still gauge-invarianf^. On fhe ofher 
hand, fhe compacf QED action infroduces phofon-phofon inferacfions which are a pure discretiza¬ 
tion effecl. However, fhere is in principle no concepfual problem in using compacf QED. In bofh 
cases (allhough if is mandatory for fhe non-compacf aclion) gauge fixing needs lo be considered. 
Indeed, wifh elecfromagnefic inferacfions one mighf be inleresled in gauge varianl amplifudes in¬ 
volving charged parficles. In fhe non-compacf case, gauge fixing is slraighlforward for Coulomb 
and Eeynman gauges [6]. Moreover, as explained in [6], fhe non-compacf acfion offers an in¬ 
teresting opporlunily for Eourier acceleration of fhe hybrid Monfe-Carlo (HMC) algorifhm. The 
argumenf goes as follows: because fhe pure gauge fheory is free, one can find a dislribulion for fhe 
HMC momenla which exaclly cancels any aufocorrelafion in fhe Markov chain. Of course Ibis is 
nol correcf anymore once quarks and gluons are coupled to fhe system. However, because of fhe 
weak coupling of QED if has been observed fhaf using Ibis parficular momenfum dislribulion slill 
considerably reduces fhe aulocorrelalions coming from EM inferacfions. 

2.3 QCD-i-QED simulations 

Up to now, essentially two approaches have been used to perform lattice QCD-i-QED simula¬ 
tions. On the one hand, one can use the so-called electro-quenched approximation which consists 
in neglecting the EM contribution to the fermionic determinant {i.e. the sea quarks are electrically 
neutral). This approach is more cost effective but it is not possible to control reliably the quenching 
effects. On the other hand, one can consider the full theory. In the past couples of year several 
groups worked or started working on such simulations. We summarize below the effort done in 
both approaches. 

2.3.1 Electro-quenched approximation 

Erom the Monte-Carlo simulation point of view, the coupling of quenched QED (qQED) to 
QCD is fairly straightforward. Eor a given lattice QCD gauge configuration, one generates a pure 
gauge QED configuration (which is simply Gaussian distributed for the non-compact lattice ac¬ 
tion). Then the QED field is used lo phase fhe QCD gauge links in fhe laffice Dirac operafor 
inverted lo obfain fhe valence quark propagalors. Anolher approach has been used in [8]: in fhaf 

^this is not the case for non-compact lattice QCD [15] 
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work the leading QED corrections are expressed as pure QCD expectation values in a perturbation 
theory fashion. In that framework, the electro-quenched approximation consist in neglecting the 
disconnected quark diagrams where the EM currents are self-contracted. 

It is easy to show that the missing contributions to the fermionic determinant are suppressed 
by both the number of colors and SU(3) flavor symmetry. Using this fact and naive dimensional 
analysis, a quantity computed in the electro-quenched approximation is expected to suffer from a 
0(10%) quenching effect relatively to its electromagnetic corrections [9]. Also, partially quenched 
chiral perturbation theory allows to provide estimations consistent with the dimensional one for the 
light meson mass splittings [16, 17]. 

A summary of the different lattice QCD-i-qQED simulation projects can be found in [18]. 
Apart from a slight update from the MIEC group [10], the actual situation is essentially identical 
concerning electro-quenched simulations. 

2.3.2 Full QCD-i-QED simulations 

There are essentially three possible ways to compute full QCD-i-QED correlation functions. 
Eirstly, it is possible to compute directly the ratio of the QCD-i-QED to QCD fermionic determinant 
for each QCD configuration. These ratios can be then used to re-weight electro-quenched data. 
This technique was first proposed in [19] and applied in exploratory calculations [20, 21]. We see 
one major limitation of re-weighting techniques applied to QCD-i-QED: as the volume gets larger, 
the computational cost of the weights increases rapidly and the signal over noise ratio decreases. 
As discussed in section 2.1, it is important to reach large physical volumes in order to control the 
large finite-volume effects generated by QED. 

A second approach already mentioned in the previous section is to perform a perturbative 
expansion in a and express the QED correction as pure QCD observables. The determinant contri¬ 
bution then appears as disconnected quark diagrams. The computation of these diagrams was never 
attempted in the context of computing EM corrections to hadronic amplitudes. However, identical 
diagrams contribute in other problems and they appear to be extremely difficult to compute. One 
can for example look at the recent study by the Mainz group [22] of the disconnected contribu¬ 
tions to the hadronic vacuum polarization. This approach has the advantage of isolating specific 
perturbative contributions which can be useful for the control and understanding of IR divergences 
in hadronic processes [23]. 

Einally, one can generate new field configurations including both QCD and QED actions in the 
HMC process. Both QCDSE [24, 25] and MIEC [14] have started simulations and BMWc achieved 
the first complete simulation program [6]. The specificities of each of the projects mentioned in 
this section is summarized in Table 2. 

3. Isospin corrections to the hadron spectrum 

In this section we review the different results concerning the isospin breaking corrections to the 
hadron spectrum. We first discuss the ambiguity of the separation of strong and EM contributions. 
Then we present the results concerning Dashen’s theorem and quark masses and finally the hadron 
mass splittings. 
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collaboration 

references 

RBC-UKQCD 

[20] 

PACS-CS 

[21] 

QCDSE-UKQCD 
[24, 25] 

BMWc 

[6] 

fermion action 

domain wall 

Wilson 

Wilson 

Wilson 

Nf 

2+1 

1 + 1 + 1 

1 + 1 + 1 

l+l+l+l 

mefhod 

re-weighfing 

re-weighfing 

HMC 

HMC 

min(M;j) (MeV) 

420 

135 

250 

195 

a (fm) 

0.11 

0.09 

0.08 

0.06-0.10 

Na 

1 

1 

1 

4 

L (fm) 

1.8 

2.9 

1.9-2.6 

CO 

1 

(N 

A^vol. 

1 

1 

2 

11 


Table 2: Summary of full lattice QCD+QED simulation programs. The MILC program [14] is too prelim¬ 
inary to know its specifications and was not included in this table. The first line is the fermion action used. 
The second line is the number of flavors used in the gauge configuration generation. The third line gives the 
simulation method used. The forth line indicates the minimum pion mass reached. The fifth line is the range 
of lattice spacing used and the sixth line indicates their number. Similarly, the seventh and eighth lines are 
respectively the range and the number of lattice spatial extents used. 

3.1 Separation of QCD and QED contributions 

In all present work, only the leading isospin corrections to hadron masses are considered. 
These corrections can be written as follows: 

Mix = aAx + SmBx + 0{a^ ,a8m,5n?) (3.1) 

where Mix is a given isospin mass splitting and 5m = mu — m^. Then it is tempting to simply 
define the leading-order QED and QCD parts of the splitting: 

Aqed^^x = o:Ax and = 5mBx (3.2) 

However, on has to be careful because of the following ambiguity: a and 5m depend on each other 
through radiative corrections. Moreover, and m^ individually depend on a with a different co¬ 
efficient because of the difference of their electric charges. Therefore, this ambiguity cannot be 
directly absorbed in higher-order isospin corrections. To make properly this QCD/QED separa¬ 
tion, one has to provide a prescription that defines fhe 5m = 0 poinf. The difference befween fwo 
prescripfions will be 0{mudOi,mud5m) up fo higher-order isospin corrections. Thus, wifh physical 
quark masses where 5m ~ Ibis discrepancy can be considered as higher-order isospin correc¬ 
tions. So as soon as a resulf is produced af fhe physical value of mud and 5m, if is reasonable fo 
consider fhaf fhe separation (3.1) is effectively unambiguous up fo higher-order 0(1%) corrections. 

Several prescripfions fo define fhe 5m = 0 poinf have been proposed in previous works. The 
mosf concepfually sfraighfforward scheme is fo renormalize fhe lighf quark masses in a given 
scheme {e.g. MS or RI-MOM) af a given scale and fo consisfenfly express every quantify as a 
function of fhese renormalized masses. This prescription was used in [26, 12, 8]. Because renor¬ 
malized quark masses can be difficull fo compute, if is also inferesfing fo consider prescripfions 
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based on hadron masses. In [9], 5m was replaced by the mass squared difference between the 
connected Tiu and dd mesons AM^ = Is it possible to show in partially-quenched chiral 

perturbation theory [16] that for physical quark masses, AM^ is directly proportional to 8m up to 
the 0(1%) higher-order corrections. In the later work [6] from the same collaboration, the pre¬ 
scription AqedAIi: = 0 was used, i.e. the £+ — mass splitting is assumed to be proportional to 
5m. If these particles would be point-like, one would have AgED^fe = 0 exactly. The authors of 
[6] found AgEoAfi: statistically consistent with 0 and no more than 0.2 MeV with AM^ = 0. 

3.2 Dashen’s theorem and light quark masses 

Dashen’s theorem corrections and individual up and down quark masses have been computed 
reliably only in the electro-quenched approximation. All existing results on Dashen’s theorem 
correction £ defined in (1.4) are presented in Figure 2. Two interesting comments can be make 
regarding these results. Firstly, although lattice results are still dominated by systematic uncertain¬ 
ties and suffer from an uncontrolled electro-quenching error, they look consistently spread around 
a common value. This contrasts significantly with the 1990’s phenomenological determinations 
of this quantity. Secondly, the lattice results seems to favor a rather large 0(70%) violation of 
Dashen’s theorem. 



Maltman and Kotchan [27] (1990) 
Donoghue ef a/. [28] (1993) 

Bijnens [29] (1993) 

Baur and Urech [30] (1996) 

Bijnens and Prades [31] (1997) 
Donoghue and Perez [32] (1997) 

Gao eta/. [33] (1997) 

Moussallam [34] (1997) 

Duncan et al. [7] (1996) (EQ, quenched) 
RBC-UKQCD [26] (2007) (EQ) 
RBC-UKQCD [12] (2010) (EQ) 

RM123 [8] (2013) (EQ) 

BMWc [35] (2014) (EQ, preliminary) 
MIEC [10] (EQ, preliminary) 


0 0.5 1 1.5 2 

£ 


Figure 2: Summary of the determination of Dashen’s theorem violation e dehned in (1.4). The green points 
represents analytical calculations from effective theory and the blue points results from lattice simulations. 
For the lattice results, the blue error bar is statistical and the red is systematic. “EQ” stands for electro- 
quenched. Results are presented in chronological order. 


Regarding the light quark mass ratio mu/ma, the lattice determinations of this number are sum¬ 
marized in Figure 3. Although there is a slight tension between the two most recent results [10, 35], 
they both agree nicely with the values from the PDG [1] and FLAG [2] reviews. The only full 
QCD-i-QED result from PACS-CS [21] seems to deviate significantly from other determinations. 
This number is the result of an exploratory calculation using re-weighting techniques and has un¬ 
known systematic errors. We believe that this effect is more a systematic effect rather than an 
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indication of a large see quark EM contribution. It is interesting to notice that all these results 
exclude strongly the = 0 solution to the strong CP problem. 


RBC-UKQCD [12] (2010) (EQ) 
PACS-CS [21] (2012) 

RM123 [8] (2013) (EQ) 

BMWc [35] (2014) (EQ, preliminary) 
MIEC [10] (EQ, preliminary) 

0.4 0.5 0.6 0.7 

mu/nid 

Figure 3: Summary of the calculations of m^lmd light quark mass ratio. Plotting conventions are identical 
to the one used in Figure 2. 



3.3 Isospin mass splittings in the hadron spectrum 

The main novelty concerning the calculation of the isospin correction to the hadron spectrum 
is the high-precision determination of isospin mass splittings of the octet baryons, the D meson 
and the newly discovered Sec. from the BMWc group [6]. These splittings where computed using 
full QCD-i-QED simulations including an active charm quark in the sea. A summary of these 
results can be found in Eigure 4. The splittings obtained in this work are in very good agreement 
with experimental values. It is interesting to notice that the S baryon splitting is obtained with a 
precision higher than the experimental measurement. Moreover, the unknown Sec splitting needed 
by charm spectrum experiments^ is predicted accurately. 

The QCDSE-UKQCD collaboration also aim at studying isospin corrections to the octet baryon 
spectrum. This group have started generating full QCD-i-QED gauge ensembles [24, 25] to deter¬ 
mine the corrections to the spectrum. The analysis is performed using the same technique based 
on SU(3) flavor symmetry as used in their previous pure QCD work [36]. The same collaboration 
also achieved the first lattice determination of the mixing between Y? and A*’ baryons [37] which 
is authorized once isospin symmetry is broken. They obtained the mass splitting between the two 
particles with a precision of 0(10%), in good agreement with the experimental value. 

Einally we summarize the theoretical determination of the nucleon splitting, including the 
QCD and QED separation, in Eigure 5. The most recent full QCD-i-QED results [6, 25] are in very 
good agreement and indicate that the crucial value of the nucleon mass splitting is indeed the result 
of a subtle cancellation between the QCD and QED contributions. 

^cf. e.g. http://www.ectstar.eu/sites/www.ectstar.eu/files/talks/after-jurgen-febl3.pdf 
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Figure 4: Summary of the results from [6]. Acg = AM^^ — AMj; + AMe is the correction to the Coleman- 
Glashow relation. 

4. Conclusion and perspective 

Lattice QCD simulations in the isospin limit are reaching a precision of 0(1%) and below on 
important observables and isospin breaking effects are becoming the dominant source of system¬ 
atic uncertainty. Therefore, it is becoming crucial to introduce isospin breaking effects in lattice 
simulation in order to provide more stringent theoretical constraints on the Standard Model. The 
main challenge in this task is the inclusion of EM interactions. 

The difficulty with QED comes from the IR singular structure of the theory. Defining correctly 
the theory in a box is non-trivial for several reasons. Eirstly, momentum quantization can introduce 
hard singularities coming from the photon field zero-mode if one uses periodic boundary condi¬ 
tions. Subtracting the zero-mode of the field is a possible solution. However, fixing the zero-mode 
constitutes a non-local constraint on the theory which can break some important properties such 
as reflection positivity. It is shown in [6] that the zero-mode subtraction proposed by Hayakawa 
& Uno [11] has a correct quantum mechanical interpretation. Beyond the zero-mode subtraction, 
QED in a finite volume suffers from large, power-like finite-size effects because of the long-range 
of the interaction. These effects are now well understood for hadron masses. The two first orders 
in the infinite-volume expansion are entirely determined by gauge invariance and can be computed 
analytically [6, 5]. Several effective theory descriptions of the higher-order, structure dependent 
corrections have been worked out [11, 12, 5, 13]. 

In its recent work, the BMWc group [6] achieved the first complete lattice calculation featur¬ 
ing full QCD and QED interactions, an active sea charm quark and a full control over the different 
sources of uncertainty. The corrections to the baryon octet and charm spectrum are obtained pre¬ 
cisely, in good agreement with experimental measurements. These simulations represent an impor¬ 
tant step toward fully non-perturbative, high precision, predictions of Standard Model observables. 
The corrections to the light meson spectrum, necessary to determine the individual up and down 
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Figure 5: Review of the theoretical determinations of the nucleon mass splitting. Results represented by a 
band correspond to works where only the QCD or QED contribution has been determined. Where possible, 
error bars are represent statistical and systematic uncertainties. The “no j3-decay” region is defined by 
M„ — Mp < me where me is the electron mass. “EQ” stands for electro-quenched. 


light quark masses are still known only through the electro-quenched approximation. 

It is now important to consider more complex quantities than hadron masses. Adding isospin 
breaking effects to the determination of hadronic decay widths is crucial to obtain high precision 
constraints on the flavor structure of the Standard Model (e.g. through sub-percent determinations 
of the CKM matrix coefficients). With EM interactions, matrix elements are significantly harder 
to determine than energy levels. Indeed, such quantity can feature IR divergences that are phys¬ 
ically cancelled by the addition of real soft photons in the final state. Recently, a proposal has 
been made [23] to deal with such divergences in the case of meson decays. Also, QCD-i-QED 
simulations can be used to perform non-perturbative computation of the hadronic corrections to 
the muon anomalous magnetic moment [42]. This quantity features an interesting discrepancy be¬ 
tween theory and experiment and it is important to reinforce the theoretical prediction to support 
the experimental effort. 
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